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Abstract 

High precision cross-section data of the deuteron-proton breakup reaction at 130 
MeV deuteron energy are compared with the theoretical predictions obtained with 
a coupled-channel extension of the CD Bonn potential with virtual A-isobar exci- 
tation, without and with inclusion of the long-range Coulomb force. The Coulomb 
effect is studied on the basis of the cross-section data set, extended in this work 
to about 1500 data points by including breakup geometries characterized by small 
polar angles of the two protons. The experimental data clearly prefer predictions 
obtained with the Coulomb interaction included. The strongest effects are observed 
in regions in which the relative energy of the two protons is the smallest. 
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The successes of meson exchange theories in the description of two-nucleon 
observables directed the research towards systems composed of three nucle- 
ons (3N), where nucleon-nucleon (NN) interaction models can be tested in 
a non-trivial environment and, moreover, additional dynamics related to the 
presence of the third nucleon can be investigated. Breakup process, with its 
variety of final states and rich information contained in the observables of 3N 
continuum, is ideally suited for such tests. Although studies of this process are 
very challenging, both in the sense of theoretical calculations as well as of pre- 
cise measurements, an important progress has taken place in this field over last 
years. At present, breakup observables can be predicted rigorously via exact 
solutions of the Faddeev equations with realistic NN potentials, combined with 
model 3N forces [1] or with the two- and three-nucleon interactions obtained by 
an explicit treatment of the A-isobar excitation [2,3,4,5] within the coupled- 
channel framework. Alternatively, the dynamics is generated by the chiral 
perturbation theory approach at the next-to-next-to-leading order [6,7,8,9,10] 
with all relevant NN and 3N contributions taken into account. In parallel, the 
data base, rather poor until recently in the region of medium energies, has 
been significantly enriched by our measurement of the ^H(d,pp)n reaction at 
the beam energy of 130 MeV. Precision of the obtained experimental data 
and coverage of a large fraction of the phase space for cross sections [11,12] 
and for analyzing powers [13,14] allows to reliably test predictions of various 
theoretical approaches. In this paper we supplement the cross-section results 
with additional configurations, characterized by small polar angles of the two 
breakup protons. 

The theoretical achievements in the 3N system have been shadowed (but also 
stimulated) by a persistent discussion about a possible bias of the obtained 
conclusions indiiccd by neglecting potentially important pieces of the inter- 
action dynamics. In all the above mentioned approaches the calculations are 
performed in a non-relativistic regime. Moreover, they ignore the long-range 
Coulomb interaction, thus can describe properly the neutron-deuteron (nd) 
system, while the data used in their verifications are precise (and numerous) 
enough only for the proton-deuteron (pd) system. Very recently the situation 
started to change. An important progress took place in both, the fully relativis- 
tic treatment of the breakup process as well as with respect to the inclusion of 
the Coulomb interaction. It is worthwhile to comment that - contrary to the 
long-trusted expectations - in both aspects the conclusions from the sector 
of elastic scattering proved not to be directly applicable for the breakup pro- 
cess, for the whole variety of its kinematical configurations. Pioneering study 
on incorporating relativity in the calculation of the elastic scattering observ- 
ables [15] indicated that those effects are indeed expected to be very small. 
However, the very recent study on implementing boost and relativistic kine- 
matics to the breakup process [16] revealed quite large (several percent) effects 
for some specific geometries. When regarding the long-range electromagnetic 
force influence, again the calculations for the elastic scattering cross section at 
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65 McV [17,18] showed an essentially negligible difference between nd and pd 
predictions, even in the cross-section minimum, the most sensitive region to 
study the 3N force effects (significant Coulomb effects are visible only at very 
small scattering angles). However, the very first set of calculations, in which 
the Coulomb effect in the breakup reaction is taken into account [19,20], indi- 
cates that they can lead to a dramatic change of the cross section magnitude 
in certain regions of phase space. In this paper the influence of the Coulomb 
interaction on the differential cross sections of the breakup reaction is further 
studied, using the extensive data set at 130 MeV deuteron energy. The aim is 
a quantitative check of the theoretical predictions which include the effect of 
the Coulomb interaction and the analysis of the dependence of this effect on 
kinematical variables. 

The theoretical predictions are based on a realistic coupled-channel poten- 
tial CD Bonn -|- A [5], allowing for a single virtual A- isobar excitation and 
thereby yielding an effective 3N force consistent with the NN force, and in- 
cluding exchanges of vr, p, oo, and a mesons. The Coulomb interaction be- 
tween charged baryons is fully included using screening and the renormaliza- 
tion approach [20,21]. The special choice of the screened Coulomb potential 
wr — w e"^''/^)", with n — A being optimal, approximates well the true 
Coulomb interaction w for distances r smaller than the screening radius R. 
Simultaneously, this potential vanishes rapidly for r > R, yielding relatively 
fast convergence of the results with respect to R and to the included partial 
waves. Because of screening, standard scattering theory is applicable and the 
three-particle transition matrices for elastic and breakup scattering, U^^\Z) 
and Uq^\z), referring to hadronic plus screened Coulomb interaction, are ob- 
tained by solving the symmetrized Alt-Grassberger-Sandhas equations [22] in 
momentum-space 



U^'^\Z)^PG^\Z) + PT^'^\Z)Go(Z)U^'^\Z) , (1) 
C/f = (1 + P)G^\Z) + (1 + P)r(^)(Z)G'o(Z)C/(^)(Z) , (2) 

using standard partial-wave basis. In Eqs. (1) and (2) Go{Z) is the free resol- 
vent, P is the sum of the two cychc permutation operators, and T^^^ [Z) is the 
two-particle transition matrix derived from nuclear plus screened Coulomb po- 
tentials; the dependence of operators on the screening radius R is notationally 
indicated. Finally, the renormalization procedure [23,24] is applied to obtain 
the scattering amplitudes in the unscreened limit. Further details are given in 
Refs. [20,21]. 

The experimental data were acquired in measurements performed at the Kern- 
fysisch Versneller Instituut (KVI), Groningen, The Netherlands. The deuteron 
beam with energy of 130 MeV was focused to a spot of approximately 2 mm 
diameter on a liquid hydrogen target of 4 mm thickness. The experimental 
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setup consisted of a three-plane multi-wire proportional chamber (MWPC) 
and of two layers of a segmented scintillator hodoscope: transmission AE and 
stopping E detectors. Position information from the MWPC was used for pre- 
cise reconstruction of the particle emission angles, while the hodoscope allowed 
to identify the particles, to determine their energies and to define trigger con- 
ditions. The AE-E wall covered a substantial fraction of the phase-space: from 
about 10° to 35° for polar angles 9 and the full (27r) range of the azimuthal 
angles 0. Registered were coincidences of the charged reaction products: the 
two protons emitted from the breakup reaction or proton and deuteron from 
the elastic scattering. More details on the experimental setup and procedures, 
as well as on the data analysis are given in Refs. [11,12]. 

In Ref. [12] high precision cross-section data of the deuteron-proton breakup 
reaction, obtained for 72 kinematically complete configurations, are presented 
for a regular grid of polar and azimuthal angles with a constant step in the 
arc-length variable S. Polar angles of the two outgoing protons, di and ^2, were 
selected between 15° and 30° with a step of 5°, and their relative azimuthal 
angle 0i2 was taken in the range from 40° to 180°, with a step of 20°. For each 
combination of the central values ^1, 62 and 0i2 the experimental data were 
integrated within the hmits of ±1° for the polar angles and of ±5° for the rela- 
tive azimuthal angle. The bin along the kinematic curve S was 4 MeV. Various 
theoretical predictions, valid for the nd system, well reproduce the experimen- 
tal dp data in several configurations. However, large discrepancies have been 
observed in geometries characterized with the smallest analyzed polar angles, 
Oi — 02 — 15°. For small 0i2 values the theories overestimate the data while 
at large 0i2 the experimental results are underestimated. Suppression of the 
pd data with respect to the nd theory at low relative azimuthal angles can be 
qualitatively understood as due to the Coulomb repulsion in the vicinity of the 
pp-FSI points. Recent results for the breakup cross sections, calculated with 
the Coulomb interaction, confirm quantitatively these expectations. Inclusion 
of the long-range electromagnetic force improves significantly the agreement 
between the theoretical predictions and the data at 9i — 62 — 15° for both, 
small and large 4>i2 values - see Fig. 1 1 of Ref. [20] . 

Sizable effects of the Coulomb interaction observed at ^1 = ^2 = 15° motivated 
us to extend the study to the lowest polar angles allowed by the detector 
acceptance, i.e. down to 12°, where the detection efficiency is still large and well 
under control. Keeping the same integration limits of ±1° for the polar angles 
and of ±5° for the relative azimuthal angle, we have analyzed 8 configurations 
with the central values of 9i — 62 — 13° and 0i2 varied from 40° to 180° 
with a step of 20°. Those configurations were not particularly interesting in 
discussing the infiuences of the 3N force since the predicted effects of the 3N 
interaction are small in this region, but they are very well suited to study 
the Coulomb effect, as can bee seen in Fig. 1. The predicted influence of 
the Coulomb force on the differential breakup cross sections is large - the 
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Fig. 1. Differential cross sections of the deuteron-proton breakup at 130 MeV 
deuteron energy, plotted as a function of the arc length S along the kinematical 
curve. The data are shown for 8 kinematical configurations characterized by the 
proton polar angles 9i = 02 = 13° and various relative azimuthal angles 4'i2, as 
indicated in the individual panels. The error bars represent statistical uncertainties 
only. Experimental data are compared to the results of calculations with the cou- 
pled-channel CD Bonn + A potential, without (dashed lines) and with (solid lines) 
inclusion of the Coulomb interaction. The results in the first panel are normalized 
to the common vertical axis by the indicated multiplier factor. 

differences between the calculations with (solid lines) and without (dashed 
lines) Coulomb force reach locally almost 25 percent. Due to high cross section 
values the statistical errors of the data points arc very small (below 3 percent). 
Discussion of the systematic uncertainties given in [12] is valid also in the 
kinematical region of the smallest polar angles, thus the systematic errors 
of the data points are on the level of 2-4 percent. Comparison of the data 
with the theoretical predictions (Fig. 1) shows that, indeed, the inclusion of 
the long-range Coulomb force in the calculations significantly improves the 
agreement with the data. Only for the configurations with equal to 80° and 
100° the calculations with the Coulomb interaction included give somewhat 
worse description of the data, while the improvement observed at the extreme 
values of 0i2 is striking. In terms of the global comparison (as discussed 
in Refs. [11,12]) of the experimental cross-section values, the inclusion of the 
Coulomb interaction for the configurations presented here (nearly 200 data 
points) leads to a decrease of the value by 38%. When all our cross-section 
data (nearly 1500 data points in 80 configurations) are used in the comparison, 
the decrease due to the inclusion of the Coulomb force amounts to 20%. 
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Fig. 2. Relative discrepancies between the experimental data and the theoretical 
predictions of the breakup cross sections as a function of the relative energy of the 
two breakup protons. The circles show the results obtained with the CD Bonn + A 
potential calculations without the Coulomb force. The full dots represent the results 
for the calculations with the inclusion of the Coulomb force. 

These results prove how important it is to include all aspects of the interaction 
dynamics in the theoretical description. 

Both, the quantitative considerations and the precise calculations indicate 
that the kinetic energy of the relative motion of the two protons, Erei, should 
be an important parameter in the studies of the Coulomb force effects. This 
quantity, equal to the total kinetic energy of the two protons in their center- 
of-mass reference frame, can be calculated on the basis of the energies and 
directions of the two emitted protons as 



rel 



\/{Ei + E. 



(3) 



where Ei,pi are respectively the total energy and momentum vector of the 
i-th {i = 1,2) proton with mass m, all expressed in the energy units. For 
each cross-section data point from the 80 kinematical configurations the cor- 
responding value of Erei was calculated and the data were sorted with respect 
to this parameter: The relative differences of the experimental and theoretical 
cross sections, {(Texp ~ crth) / {<^th) , were determined and plotted as a function of 
E'rei ill Fig. 2. The relative differences have been calculated using theoretical 
prediction without (circles) and with (full dots) the long-range electromag- 
netic force. As expected, the strongest influence of the Coulomb interaction 
can be observed at the smallest values of Ej-ei- In that region inclusion of the 
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Coulomb force strongly improves the agreement between the data and the the- 
oretical description, though the discrepancies are not completely removed. In 
the medium Erei range, 8-18 MeV, a perfect consistency between the data and 
the results of the calculations with the Coulomb force included is reached. At 
higher relative energies the theoretical calculations tend to underestimate the 
data, though it is worthwhile to note that for Erei above 20 MeV the influence 
of the Coulomb force becomes practically negligible. This observation can be 
used to argue that the nd calculations can be safely applied to the pd breakup 
data in the regions of sufficiently large Ej-ei values. 

Theoretical calculations using the coupled-channel CD Bonn + A potential 
predict dramatic changes of the breakup cross section distributions due to 
Coulomb in configurations characterized by very small polar and relative az- 
imuthal angles between the two outgoing protons. This statement is partially 
supported already by inspecting the first panel in Fig. 1. To demonstrate a still 
stronger action of the Coulomb force we analyzed in addition a kinematical 
configuration lying at the very edge of the experimental acceptance, with the 
central angular values of 9i = 02 = 13° and (pu = 20°. The event integration 
ranges have been kept as in the previous cross section analysis. It should be 
pointed out that at so small relative azimuthal angle between the two pro- 
tons the detector acceptance is significantly reduced due to its granularity. 
Events in which both breakup protons hit the same detector (AE or E) are 
lost due to uncertain energy information. The cross-section results have to be 
corrected for these losses. The correction factors are determined on the basis 
of the GEANT simulation (more details can be found in [12]), however, when 
the effects are large, the acceptance losses are sensitive to even small geomet- 
rical inaccuracies of the setup, what can affect the cross section normalization. 
Therefore, the overall systematic uncertainty of the cross-section data for the 
here discussed configuration is estimated at the level of about 6-8 percent. The 
result is presented in Fig. 3. In this spectacular case the Coulomb interaction 
leads not only to a strong suppression of the cross section (predictions disre- 
garding the Coulomb force reach a value of 12 mb-sr~^-MeV^^), but also to a 
distortion of the distribution, causing a dip at the minimal relative energy of 
the two brakup protons. This behavior is very well seen in the data as well. 

All the performed comparisons, local and global, show that our cross-section 
data are much better desribed by the predictions in which the long-range 
Coulomb force is taken into account. It can be also seen that across the breakup 
phase space there are regions in which the sensitivity to the Coulomb force is 
practically negligible. In these regions one may use nd calculations to compare 
with pd data, but elswhere the pd calculations with the full treatment of 
the long-range Coulomb force are paramount to the proper interpretation of 
the experimental pd results and conclusions on the underlying force models. 
Inclusion of the Coulomb interaction allows now a cleaner theoretical insight 
at possible shortcommings of the models of the hadronic dynamics. 
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Fig. 3. Differential cross section of the breakup reaction for the kinematical con- 
figuration characterized by the proton polar angles 9i = 02 = 13° and their relative 
azimuthal angle ^12 = 20° (dots) . The solid line represents the coupled-channels pre- 
dictions of the CD Bonn -|- A potential including the Coulomb force. The dashed 
line shows the results of calculations disregarding the Coulomb force - truncated 
vertically for a better insight into the cross section behavior in the central S region. 
The dash-dotted line and the righ-hand side scale present the dependence of the 
relative energy of the two breakup protons along the S-axis. 



This study makes an important step towards a precise and complete descrip- 
tion of the breakup observables, which should eventually include all aspects 
of the medium-energy reaction mechanism. The theoretical predictions show 
that the effects of the Coulomb force, relativity and of the 3N interaction 
affect the breakup observables in different ways and with varying strength 
when inspecting the full reaction phase space. Such selectivity makes possible 
tracing the details of certain effects in regions where the others are proved to 
have relatively small inffuences. This is e.g. true for studies of the 3N forces - 
even if the Coulomb effects in the breakup cross sections are large, there are 
regions in which their influence is much smaller than the expected effects of 
the additional nuclear dynamics. With even larger experimental coverage of 
the breakup phase space with respect to several observables and for various 
beam energies, the eventually established pattern of discrepancies between the 
data and the calculations might help to improve the understanding of the full 
dynamics of the 3N system. 
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